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We present a detailed investigation of the magnetic and superconducting properties of 
Cai_ 2 ;Na 2 ;Fe 2 As 2 single crystals with x = 0.00, 0.35, 0.50, and 0.67 by means of the local probe 
techniques Mossbauer spectroscopy and muon spin relaxation experiments. With increasing Na- 
substitution level, the magnetic order parameter as well as the magneto-structural phase transition 
are suppressed. For x = 0.50 we hnd a microscopic coexistence of magnetic and superconducting 
phases accompanied by a reduction of the magnetic order parameter below the superconducting 
transition temperature Tc. A systematic comparison with other 122 pnictides reveals a square-root 
correlation between the reduction of the magnetic order parameter and the ratio of the transi¬ 
tion temperatures, T^/Tn, which can be understood in the framework of a Landau theory. In the 
optimally doped sample with Tc « 34 K, diluted magnetism is found and the temperature de¬ 
pendence of the penetration depth and superfluid density are obtained, proving the presence of two 
superconducting s-wave gaps. 

PACS numbers: 74.70.Xa, 76.75.-|-i, 76.80.-|-y, 74.62.Dh 


I. INTRODUCTION 

Since the discovery of superconductivity in iron 
pnictidestheir electronic phase diagrams, characterized 
by a close proximity of magnetic and superconducting 
phases, have been explored in great detail. Of particular 
interest are the regions of the phase diagrams showing 
a crossover from magnetic order to superconductivity. 
Magnetic spin fluctuations, enhanced in the vicinity of 
a magnetic quantum critical point, can play an impor¬ 
tant role in the formation of Cooper pairs. In addition, 
the interplay of both magnetic order and superconduc¬ 
tivity can lead to a phase with microscopic coexistence 
of both ground statesi^^— As both states compete for the 
same electrons at the Fermi surface, the magnetic order 
parameter may be reduced below TcA— In this work, we 
studied the system Cai_a,Naa;Fe 2 As 2 . The parent com¬ 
pound CaFe 2 As 2 shows spin density wave order below the 
Neel-temperature = 165 The magnetic phase 
transition is accompanied by a structural phase transition 
from a tetragonal to an orthorhombic structurei^ Increas¬ 
ing the Na amount in Cai_a,Naa;Fe 2 As 2 , the magneto- 
structural phase transition is suppressed until it vanishes 
at a critical Na concentration x ss 0.35ji Superconduc¬ 
tivity is found for x >0.3 with a maximum of ~ 34 K 
at an optimal doping of x ~ 0.66. However, the interac¬ 
tion of superconductivity with the magnetic spin density 
wave in the region of 0.3 < x < 0.35 has not been 

conclusively determined. 

We studied the magneto-structural phase transition 
as well as both the superconducting and the mag¬ 
netic order parameter and their interaction for differ¬ 
ent Na-substitution levels. We find a suppression of the 
magneto-structural phase transition upon Na substitu¬ 


tion. For X = 0.35 and 0.50, microscopic coexistence of 
magnetic order and superconductivity is observed. For 
the latter Na-substitution level, a reduction of the mag¬ 
netic order parameter is observed below the supercon¬ 
ducting transition temperature. For x = 0.67, two super¬ 
conducting gaps with s-wave symmetry can be deduced 
from the temperature dependence of the superfluid den¬ 
sity. 


II. EXPERIMENTAL 

We examined mosaics of Cai_xNaa;Fe 2 As 2 single crys¬ 
tals, which were grown by the self-flux technique as de¬ 
scribed by Johnsten et al^ The samples were character¬ 
ized by energy-dispersive X-ray spectroscopy (EDX), X- 
ray diffraction (XRD), susceptibility, magnetization, and 
specihc-heat measurements. The stoichiometry of the 
examined samples is a; = 0.00, 0.35, 0.50, and 0.67 as 
determined by EDX. A characterization of the magnetic 
properties was performed using SQUID magnetometry in 
large and small external magnetic fields. Mossbauer spec¬ 
troscopy (MBS) experiments were performed in transmis¬ 
sion geometry in a temperature range between 4.2 and 
300 K using a CryoVac Konti IT cryostat. As the 7 
source, a ®^Co in rhodium matrix was used. The sin¬ 
gle crystals were aligned with the crystallographic c-axis 
along the 7 direction. To analyze the data, the hyperfine 
Hamiltonian including electric quadrupole and magnetic 
hyperfine interactions was diagonalized. In the paramag¬ 
netic temperature regime, the spectra were described by 
a doublet pattern, whereas in the magnetically ordered 
state a sextet pattern was used. The magnetic order pa¬ 
rameter is deduced from the ®^Fe magnetic hyperhne held 
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B\ii- The isomer shift 5 is given with respect to a-Fe. 

Muon spin relaxation (^SR) experiments were per¬ 
formed at the ttMS and PiEl beamlines of the Swiss 
Muon Source at the Paul Scherrer Institut, Switzerland, 
using the GPS and DOLLY spectrometers. The single 
crystals were aligned with the crystallographic c-axis 
along the muon beam. Positively charged muons 
which are nearly 100 % spin-polarized due to parity vio¬ 
lation during the pion decay, are implanted in the sample 
and thermalize at interstitial lattice sites, where they ra- 
dioactively decay with a lifetime of 2.2 ^s into two neu¬ 
trinos and one positron. As the muon decay involves the 
weak interaction, where parity conservation is violated, 
the positron is predominately emitted along the direction 
of the muon spin at the moment of the decay. Measur¬ 
ing the time-resolved angular distribution of the emit¬ 
ted positrons allows to extract the time evolution of the 
muon spin polarization P{t). The initial muon spin was 
rotated by approximately —45° with respect to the beam, 
which allows to measure the time evolution of the muon 
spin polarization T c and || c by analysing the asymme¬ 
try of the up-down and forward-backward detector pair 
respectively. If not stated otherwise, all presented mea¬ 
surements refer to the up-down detector pair. The muon 
spin relaxation was measured for temperatures ranging 
from 1.6 K up to 300 K in zero field (ZF) and transverse 
magnetic fields (TF) up to 130 mT. The /xSR data were 
analyzed using the MUSRFIT software packagei^^ 

In a magnetically ordered material, the muon spin ex¬ 
hibits a Larmor precession with a frequency which 
is related to the local magnetic field B at the muon 
site by = B^^/{2tt) (muon gyromagnetic ratio 7 ^ = 
27r X 135.5 MHz/T). The muon spin precession can be 
described in single crystals using the functionii 
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with N denoting the number of inequivalent muon sites 
contributing to the ^SR signal, weighted by P{vi). In 
the case of 100 % magnetic ordering, = 1 - 

describes the oscillating part of the signal. In contrast 
to powder samples, where Ai = 2/3 due to spatial aver¬ 
aging, for single crystals Ai G [ 0 , 1 ]. At, describing the 
damping of the oscillation, is a measure of the width of 
the static field distribution, also including dynamic con¬ 
tributions due to magnetic fluctuations. The damping of 
the non-oscillating part, described by At, is caused by 
dynamic magnetic fluctuations only»i^ 

To study the superconducting properties of the sam¬ 
ple with //SR, a magnetic field /toLlext was applied 
parallel or perpendicular to the muon beam. In 
type-II superconductors, a vortex lattice is formed for 
fJ-oHci < ^J'oHext < ^J'oHc 2 resulting in a spatial mag¬ 
netic field distribution^^ This magnetic field distribu¬ 
tion causes an additional damping of the muon spin os¬ 


cillation, which can be modeled via additional Gaussian 
terms of the formic 
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P{t) = ^ ^P^ cos{'jf^BA + 
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where cr/ describes the damping due to superconductivity, 
CTjv the damping due to nuclear magnetic dipol contribu¬ 
tions and ip denotes the angle between the initial muon 
spin direction and the positron detector. 

The second moment (AR^) of the internal magnetic field 
distribution, n(B), is given bjd^ 
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and related to the London penetration depths A by the 
relationi^ 


(AR2)= 0.00371 $ 2 ^, (4) 

where $0 denotes the magnetic flux quantum. A is related 
to the Gooper pair density by rig oc 1/A^. 

Therefore, /tSR allows to independently measure the 
magnetic and superconducting order parameters via the 
determination of the zero field muon spin precession fre¬ 
quency Vi and the London penetration depth, respec¬ 
tively. Moreover, //SR is able to distinguish between non¬ 
magnetic and magnetic volume fractions in the sample. 


III. RESULTS AND DISCUSSION 

Measurements of the magnetic susceptibility of 
GaFe 2 As 2 are reported by Harnagea et a/.— Magnetic 
susceptibility measurements in an applied field of 1 T 
parallel to the a 6 -plane are shown in Fig.jT] For x = 0.35, 
0.50, and 0.67 a nearly linear decrease of the magnetic 
susceptibility is observed in the paramagnetic region, 
which is found also for many other iron pnictidesj^^— 
The kink at 143 K, 119 K, and 150 K for t = 0.35, 0.50, 
and 0.67, respectively, indicates the onset of the anti¬ 
ferromagnetic (AFM) ordering. Magnetic susceptibility 
measurements in an applied field of 2 mT parallel to the 
afe-plane are shown in Fig.|2| For x = 0.35, below 17 K a 
broad superconducting transition occurs. As the diamag¬ 
netic shielding is not fully developed, only parts of the 
sample show superconductivity. For x = 0.50, a two-step 
superconducting transition occurs. Below 34 K, a slightly 
negative magnetic susceptibility is measured, which indi¬ 
cates a superconducting phase in a small volume of the 
sample. Below 17 K, a broad second transition occurs, 
where bulk superconductivity is formed resulting in the 
full superconducting response. For the further treatment 
of the X = 0.50 sample, 17 K will be considered as the 
superconducting transition temperature. This two-step 
behavior as well as the broad transition indicates an in- 
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FIG. 1. Temperature dependence of the static susceptibility 
X = M/H of Cai_a;Naa;Fe 2 As 2 for x = 0.35, 0.50, and 0.67. 
The measurements were performed at an applied field of 1 T 
parallel to the afe-plane. denotes the magnetic transition 
temperature. 

homogeneous sample. The sample with x = 0.67 shows 
bulk superconductivity below 34 K. The sharp transition 
indicates a homogeneous sample. 
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FIG. 2. Temperature dependence of the magnetic susceptibil¬ 
ity Xv [zero (ZFC) and field cooled (FC)] of Cai_ 2 ;Naa;Fe 2 As 2 
for X = 0.35, 0.50, and 0.67. The measurements were per¬ 
formed at an applied field of 2 mT parallel to the a6-plane. 


A. The parent compound CaFe2As2 

The time evolution of the muon spin polarization in 
ZF is shown in Fig. [3] In the paramagnetic temper¬ 
ature regime, a weak Gauss-Kubo-Toyabe damping of 
the signal is observed caused by the dipole-dipole in¬ 
teraction of the muon magnetic moment with randomly 
oriented nuclear magnetic moments. The temperature 
dependence of the magnetic volume fraction (MVF) is 
shown in Fig. 0] To describe MVF as a function of 
T, two temperatures are defined: T™®®* describes the 
highest temperature with a finite MVF and de¬ 

scribes the highest temperatures with MVF = 100 % = 1. 
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FIG. 3. ZF-pSR time spectra of Gai_a,Naa;Fe 2 As 2 for charac¬ 
teristic temperatures in the paramagnetic and magnetically 
ordered state. The well-defined muon spin precession for 
X = 0.00, 0.35, and 0.50 indicate long-range commensurate 
magnetic order. For x = 0.67, non-magnetic behavior is found 
down to lowest measured temperatures. 

The sharp transition between = 167 K and 

jimo % _ indicates a homogeneous sample. Below 

167 K, two magnetically inequivalent muon stopping sites 
A and B with a temperature-independent occupation ra¬ 
tio of Pa-Pb= 80:20 were observed as it was found in 
BaFe 2 As 2 »^ The signal fraction corresponding to muons 
stopping at site A show a well-defined sinusoidal oscilla¬ 
tion below 167 K. This indicates static long-range com¬ 
mensurate magnetic order. The temperature dependence 
of VA is shown in Fig. [5] The step-like behavior indicate 
a first-order transition, as it was seen in SrFe 2 As 2 i^ The 
signal fraction corresponding to muons stopping at site 
B show an exponential relaxation below 167 K and a 
well-defined sinusoidal oscillation below 60 K. This in¬ 
dicates a a broad held distribution at temperatures be¬ 
tween 60 K and 167 K at the muon stopping site B sup¬ 
pressing a coherent oscillation of the muon spins, which 
is contrary to the observations in BaFe 2 As 2 , where the 
two oscillation frequencies were obtained at all tempera¬ 
tures below Tat,— but consistent with LaOFeAs.— Below 
60 K, the two precession frequencies have a temperature- 
independent ratio of vaI^b ~ 1-9. Therefore, the mag¬ 
netic held at the muon stopping site A is higher than on 
site B. Site A is located next to the FeAs-layer.— The 
smaller value of the magnetic held at site B indicates 
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FIG. 4. Magnetic volume fraction as a function of temper¬ 
ature (lines are guides to the eye only) of Cai_ 2 ;Naa:Fe 2 As 2 
obtained from ZF-^SR and Mossbauer spectroscopy for x = 
0.00, 0.35, and 0.50. The inset shows the temperature depen¬ 
dence of the Mossbauer linewidth w. The abrupt increase of 
w indicates magnetic ordering. 


a muon stopping site more distant from the FeAs layer. 
As the precession frequency is proportional to the mag¬ 
netic field at the muon site, this ratio is different from 
BaFe 2 As 2 ^ and SrFe 2 As 2 )^ showing ratios of 4.1 and 
3.4. The ionic radii of the alkaline earth metals scale 
like Ca < Sr < Bai^ As a consequence, the crystal¬ 
lographic c-axis is shortest for CaFe 2 As 2 i and longest 
for BaFe 2 As 2 i^ By shrinking the crystallographic c-axis, 
the distance between site B and the ordered iron mag¬ 
netic moments in the FeAs-layer is reduced. The muon 
spin interacts with the ordered electronic moments via 
dipole-dipole and transferred Fermi-contact interaction. 
Both interactions are sensitive to the distance between 
the muon spin and the iron ordered moments. This im¬ 
plies that the change of the frequency ratio in the un¬ 
doped compounds has a structural origin. 

Mossbauer measurements by Alzamora et al^ on 
CaFe 2 As 2 show a first-order-like magneto-structural 
phase transition below 173 K. They report a saturated 
magnetic hyperfine field of ss 10.1 T at 4.2 K. The angle 
between the magnetic hyperfine field and the principal 
axis of the electric field gradient was reported as 94(4)°. 
Therefore, as the principal axis of the electric field gradi¬ 
ent is parallel to the crystallographic c-axiS)^ the mag¬ 
netic hyperfine field is located in the a6-planei^i21 These 
Mossbauer results are consistent with our /iSR results 
showing that the magnetic moments are located in the 
a6-plane. 


B. Underdoped Cao.65Nao.35Fe2As2 and 
Cao.5oNao.5oFe2As2 

Mossbauer spectra for characteristic temperatures in 
the paramagnetic and magnetically ordered states are 
shown in Fig. [51 In the paramagnetic state, an asym¬ 


metric doublet structure, which is caused by the interac¬ 
tion of the nucleus with an electric field gradient (EFG), 
was observed for both stoichiometries. However, in the 
principal-axis system, the EFG is fully determined by 
its 2 -component Vzz and the asymmetry parameter rj 
(which turned out to be zero for all investigated tem¬ 
peratures). At room temperature, for x = 0.35 and 
0.50 a value of = 11.2(5) V/A^ was obtained. 

With decreasing temperature, ^22 slightly increases to 
13.0(5) V/A^ directly above the magnetic transition tem¬ 
perature as shown in Fig. [T] With the onset of magnetic 
order, Vzz increases to 19(2) V/A^ and 18(3) V/A^ for 
X = 0.35 and 0.50, respectively, and remains constants 
within error bars down to lowest temperature. This in¬ 
crease can be assigned to a magneto-structural phase 
transition as observed for other 122 coinpoundsj ^^'^^ Fur¬ 
thermore, the principal axis of the EFG is parallel to the 
crystallographic c-axis. The temperature dependence of 
the angle 9 between the principal axis of the EFG and 
the magnetic hyperfine field, obtained from a fit to the 
data, is shown in Fig. [71 We have obtained 9 = 80(5)° 
and 71(5)° for x = 0.35 and 0.50, respectively, in the 
fully ordered state were obtained. These values indicate 
a tilting of the magnetic moments out of the a6-plane. 

The time evolution of the muon spin polarization in 
ZF on Gai_a,Naa;Fe 2 As 2 (with x = 0.00, 0.35, 0.50, and 
0.67) is shown in Fig. |31 As for the parent compound 
described above, a weak Gauss-Kubo-Toyabe damping 
of the signal is observed caused by the dipole-dipole in¬ 
teraction of the muon magnetic moment with randomly 
oriented nuclear magnetic moments in the paramagnetic 
temperature regime. The onset of long-range commen¬ 
surate magnetic ordering below = 160(2) K and 

125(3) K for 3 ; = 0.35 and 0.50, respectively, is indi¬ 
cated by the appearance of a well-defined muon spin pre¬ 
cession with two frequencies va and vb- The occupa¬ 
tion probability is independent of temperature and the 
Na-substitution level {P{va)'-P{i^b)= 80:20). We con¬ 
clude that two magnetically inequivalent muon sites are 
present in a homogeneous magnetically ordered phase. 
The temperature dependence of the muon frequency va 
for both samples is shown in Fig. jS^a) andjSKb). The 
onset temperature is consistent with the results of the 
macroscopic magnetizations measurements by SQUID 
magnetometry. The two frequencies have a temperature- 
independent ratio of vaI^b ~ 4 and 8 for a; = 0.35 and 
0.50, respectively, which are different compared to the 
undoped compound showing vaI^b ~ 1-9. This change 
is more drastic compared to other 122 compounds. 
BaFe 2 As 2 shows a ratio of 4.1;^ which increases to 4.47 
for Bao.77Ko.23Fe2As2^ and 4.5 for Bao.7Nao.3Fe2As2i^ 
SrFe 2 As 2 shows a ratio of 3A^ which changes to 3.6 for 
Sro. 5 Nao. 5 Fe 2 As 2 <^ Substituting an alkaline-earth metal 
by an alkaline metal with a smaller ionic radius (Ba^Na, 
Sr—>Na), leads to a smaller change in the c-axis param¬ 
eter and the change of the precession frequency ratio 
vaI^b occurs at lower substitution levels than by sub¬ 
stituting with an alkaline metal with a larger ionic ra- 
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FIG. 5. Temperature dependence of the magnetic order parameter for (a) x = 0.00 and 0.35 and (b) x = 0.50 including best 
order parameter fits according to Eq. The magnetic hyperfine field for x = 0.50 is Gaussian-distributed. The average value 
of the Gaussian distribution is shown with one standard deviation as error bar. The inset in (b) shows the low temperature 
regime, where superconductivity occurs below Tc = 17 K. 



FIG. 6. Mossbauer spectra of Cai_a:Naa:Fe 2 As 2 for charac¬ 
teristic temperatures in the paramagnetic and magnetically 
ordered state. The sextet structure of the Mossbauer spec¬ 
tra for X = 0.35 and 0.50 indicate long-range commensurate 
magnetic order. For x = 0.67, a broadening below « 150 K 
is observed, which cannot be described by a change in Vzz 
indicating small magnetic fields at the iron nucleus. 


dius (Ca^Na, Ba^K, Sr^K) 

scale like Ca < Na < Sr < Ba < Kj^ By increasing the 
Na-substitution level in Cai_a;Naa;Fe 2 As 2 , the crystal¬ 
lographic c-axis is elongated resulting in an increased 
distance of the muon spin at the site B from the iron or¬ 
dered moments in the FeAs-plane. In addition, the small 
tilting of the magnetic moments out of the a6-plane may 
also change the magnetic field at the muon site. How¬ 
ever, a change in the Fermi-contact interaction for both 
muon stopping sites cannot be ruled out. Additionally, 
a signal fraction in the forward-backward detector pair 
(not shown here) showing an exponential relaxation is 
observed in the magnetically ordered phase. This sig¬ 
nal fraction increases as a function of Na-substitution 
level. This is consistent with the tilting of the magnetic 
moments out of the a6-plane as observed by Mossbauer 
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FIG. 7. Temperature dependence of 0, the polar angle be¬ 
tween the magnetic hyperfine field Bm and the principal axis 
of the electric field gradient, Vzz- The increase of Vzz in¬ 
dicates a change in the electromagnetic environment of the 
®^Fe-nucleus and corresponds to the magneto-structural phase 
transition. The vertical dashed lines denotes the highest tem¬ 
perature, where the samples shows a MVF of 100 %. 


spectroscopy. 

The temperature dependence of the magnetic volume 
fraction was determined independently by ZF- and TF- 
/rSR measurements as well as Mossbauer spectroscopy. 
TF-/rSR experiments were performed by applying a mag¬ 
netic field of 5 mT perpendicular to the initial muon spin 
polarization. The temperature dependence of the MVF 
is shown in Fig. S) The onset of the magnetic ordering is 
also indicated by the abrupt increase of the Mossbauer 
line width w due to appearance of a magnetic hyperfine 
field. This increase is shown in the inset of Fig. plead¬ 
ing to Tonset = 161(2) K and 125(3) K for a; = 0.35 
and 0.50, respectively. Therefore, the obtained charac¬ 
teristic temperatures of the magnetic phase transition 
are of equal value within error bars for both /iSR and 
MBS and coincide with the magnetic phase transition 
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temperature obtained by magnetic susceptibility mea¬ 
surements, T^. The temperature width of the phase 
transition, AT = , increased from 4 K 

to 21 K and 45 K for a; = 0.00, 0.35, and 0.50, respec¬ 
tively. We attribute this increased width of the mag¬ 
netic phase transition to the increased degree of disorder 
due to the Na-substitution. The MVF is constant below 
T]^° = 163(2) K, 140(3)K and 80(3)K for a; = 0.00, 

0.35, and 0.50, respectively, proving bulk magnetic order. 
In the magnetically ordered state, a well-resolved sextet 
was observed in the Mossbauer spectra for x = 0.35, as 
it is shown in Fig. [6] at T = 4.5 K. This proves a static 
commensurate magnetic ground state with a well-defined 
hyperfine field. For x = 0.50, the sextet is less clearly re¬ 
solved and the spectra were modeled using a Gaussian 
distribution of magnetic hyperfine helds. This takes into 
account a higher degree of disorder than in the sample 
with X = 0.35. Consistently also the /iSR transverse re¬ 
laxation rate is largest for x = 0.50 as can be seen by 
a much faster suppression of the ZF oscillation in Fig. [3] 
compared to lower Na-substitution levels. The tempera¬ 
ture dependence of the obtained magnetic hyperfine field 
of X = 0.35 is shown in Fig. I5al The temperature de¬ 
pendence of the mean value of the obtained Gaussian 
distributed magnetic hyperfine field i?hf for x = 0.50 is 
shown in Fig. I5bl 

Both VA and B\a were analyzed using a fit to the 
temperature-dependent order parameter (M) of the form 


the magnetic order parameter of « 65 % appears. 
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FIG. 8. Low-temperature saturation values of the magnetic 
hyperfine field Bhf, the muon spin precession frequency va 
and 9, the angle between Bhf and the principal-axis of the 
EFG, as a function of the Na-substitution level. Upon Na 
substitution, an reduction of the magnetic hyperfine field and 
the muon spin precession frequency is observed proving a re¬ 
duction of the magnetic order parameter. The decrease 6 in¬ 
dicate a tilting of the magnetic moments out of the a6-plane 
as a function of the Na-substitution level. 
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C. Landau theory for coupled order parameters 


M{T) = M{T 



( 5 ) 


for X = 0.35 and 0.50, respectively, above the supercon¬ 
ducting transition temperatures (17 K in both cases). 
For X = 0.35, a fit to Eq.® above Tc for both i^a{T) 
and B\i{{T) represents the data in the whole temper¬ 
ature range, i.e., also below the superconducting tran¬ 
sition. Therefore, no interaction between the magnetic 
and superconducting order parameter is detectable. For 
X = 0.50, the muon frequency is reduced by approxi¬ 
mately 7 % below 17 K while the magnetic hyperfine field 
shows no reduction and is well-represented by Eq.(IS|). 
This reduction of the muon precession frequency proves 
the microscopic coexistence of magnetic order and su¬ 
perconductivity and their competition. Possibly B\^{{T) 
shows no signatures for a reduction below 17 K, since the 
reduction of 7 % is within the hyperfine field error bars. 
Alternatively, the muon precession frequency may be re¬ 
duced due to a spin reorientation below Tc rather than 
a reduction of the magnetic order parameter. However, 
in this scenario the angle 6 between the magnetic hyper¬ 
fine field Rhf and the principal axis of the EEG would 
change. The temperature dependence of 6 is shown in 
Fig.m Since there is no systematic change of 9 below 
Tc observed within error bars, a significant spin reori¬ 
entation can be ruled out. Both order parameters are 
weakly coupled compared to Bai_a;Naa;Ee 2 As 2 ,— where 
a reduction of the muon spin frequency and therefore of 


In the case of coexistence of magnetic order and su¬ 
perconductivity, an interaction between both order pa¬ 
rameters is expected to be present. This may change the 
magnitude of the order parameters and alter the crit¬ 
ical temperatures with respect to the decoupled situa¬ 
tion. We used a Landau theory to describe the coupling 
between the superconducting order parameter '0(7*) and 
the magnetic order parameter M{r) in the case of coex¬ 
istence and to describe the dependence of the reduction 
of the magnetic order parameter on the critical temper¬ 
atures. The homogeneous free-energy functional in the 
absence of an external field is given 


F[0, M] = / d?r 


+ ^l'0(F)|'^ + ^|M(r)|^ 


where the coupling between of the two order parameters 
is contained in the last term. Following the common 
approach, a and a are described as 

a =ao(T — Tato), (7) 

a =ao{T - Tco), (8) 

where Tato and TcO denote the bare magnetic and su¬ 
perconducting transition temperatures, respectively. The 
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bare transition temperatures describe the decoupled case. 
To ensure that the free-energy has a minimum and an 
that the two order parameters compete with each other, 
f3 > 0, b > 0, and d > 0 are required. The two 
order parameters are obtained by minimizing the free- 
energy functional F[ijj, M], The order parameters in the 
coexistence region are then given by 

ab-ad < 0, (9) 

= - 6^ - < 0. (10) 

This solution is stable for a sufficiently small coupling 
d, which satisfies bfd — d^ > 0. For the further treat¬ 
ment, Tco < Tnq is assumed, which is the case for 
sufficiently low doping in the 122 compounds. This re¬ 
sults in Tjvo = Tat. At Tc, where = 0, Eq. dH]) 
reduces to 


|^(Tc)P = 0 = ao6(T, - r,o) - aod(T, - T^) (11) 


which results in a linear relation between the measured 
and bare superconducting transition temperatures. The 
temperature dependence of the magnetic order parame¬ 
ter in the coexistence regime is then given by 


MI{T) = 


1 


bp-d? 


[a^d - ao/3]r 


a ,2 

“OP - 

0 


Tn- 


1 - 


ood 

aob 


Tc 


( 12 ) 


To investigate the reduction of the magnetic order pa¬ 
rameter below the superconducting transition tempera¬ 
ture, the ratio of \Mco\^{T) and \Mq\'^{T) = aQ{Tj^ — 
T)/biorT < Tc is calculated. The reduction is maximal 
for T —0. The ratio is then given by 


^(T = O') = 1 - ^ 0^0^ ~ “0^ 
m2 ao 6/3 - d^ Tn 


(13) 


Literature data for |Mcop/|Mop(T = 0) as a function 
of Tc/Tfi for various 122 compounds, measured by neu¬ 
tron scattering and fxSK, are shown in Fig. [ ^ 23 , 30 , 37- _^ 
In addition, the reduction of the structural order pa¬ 
rameter S = (a — b)/{a -I- 6) is plotted in Fig. [HI 
The magnetic and structural order parameters show 
the same behavior, as they are strongly coupled in the 
122 pnictides3iSli29£-^ |Mco|Vl^oK(7’ = 0) shows a 
nearly linear decrease as a function of T^/Ttv and hence 
the reduction of the magnetic order parameter shows 
a square root behavior. This corresponds to a con¬ 
stant slope in Eq. m, which implies a similar coupling 
strength of the magnetic and superconducting order pa¬ 
rameters in the 122 compounds. For Tc/Tm ^ 1, ex¬ 
perimental data show no reduction of the magnetic or¬ 
der parameter in the coexistence regime. This may re¬ 
sult from the fact that the samples show SC only in 


parts of the sample volume, which has been observed 
for, e.g., Cao,65Nao.35Fe2As2. Equation (US qualitatively 
describes the reduction of \Mco\^/\M q\^[T = 0) with in¬ 
creasing Tc/Tn ratio for 0 < Tc/T^ < 0.7. The system¬ 
atic deviations for Tc/Ti^ > 0.7 indicate an increase in 
the coupling strength. This means that the coupling is 
more effective if Tc ~ Tjv. 


D. Magneto-structural phase transition 

The lattice dynamics of Cai_ 2 ;Naa;Fe 2 As 2 was inves¬ 
tigated by analysing the temperature dependence of the 
Mossbauer-Lamb factor. The recoilless fraction / was ex¬ 
tracted from the Mossbauer spectra using the absorption 
area methodi^ In the Debye-approximatiou)^ 


/ oc exp 




l-b4 




dx 


(14) 

with Boltzmann constant fcs, Debye temperature 9d and 
the recoil energy Er. The temperature dependence of 
the relative recoilless fraction /(T)//(4.2 K) is shown 
in Fig. [TUI The data for temperatures above and be¬ 
low the phase transition were analyzed using Eq. ([T4l) to 
investigate the influence of the magneto-structural tran¬ 
sition on /. We accordingly define Debye temperatures 
0™ and 0^^. The step at 151 K for x = 0.35 and the 
gradual increase between 60 K and 125 K for a; = 0.50 
are attributed to an effective absorber-thickness effect. 
The magnetic phase transition leads to an increase of 
the absorption area compared to the paramagnetic re¬ 
gion, due to a splitting of the resonance lines.— The 
obtained Debye temperatures are shown in Tab. HI As the 


X 

0™/K 

0 AFM/K 

0.00 

273 M 

271^ 

0.35 

200(32) 

208(16) 

0.50 

203(12) 

219(28) 

0.67 

213(6) 



TABLE I. Debye temperatures Od obtained using Eq. (HI 
above (PM) and below (AFM) the magnetic phase transition. 

Debye temperatures 0™ and 6^^ does not change sig¬ 
nificantly within error bars at the phase transition, we 
conclude that the lattice dynamics do not change upon 
the magneto-structural phase transition. Also, the ob¬ 
tained Debye temperatures are constant within error bars 
for all X. Therefore, the lattice dynamics are independent 
of the Na-substitution level in the investigated substitu¬ 
tion range 0.35 < x < 0.67. In comparison with the 



























FIG. 9. |McoP/|A?oP(r = 0) as a function of T^/Tn for various 122 compounds.— |McoP denotes the magnetic order 
parameter in the region of coexistence with superconductivity. |Mop denotes the magnetic order parameter without any 
superconductivity. The reduction of the magnetic order parameter in the coexistence region increases with increase of T^/Tn 
following Eq. (O, the best linear fit fot T < 0.7 is shown as the solid line. 



temperature (K) 


FIG. 10. Temperature dependence of the relative recoilless 
fraction /(r)//(4.2 K) including fits following Eg. ITU 


undoped compound CaFe 2 As 2 with 9^^ r; 270 the 
lattice is softened considerably. 

To further study the magnetic properties at the 
magneto-structural phase transition, the temperature de¬ 
pendence of the isomer shift 5 was analyzed. (5 is a mea¬ 
sure for the electron density at the Fe-nucleus. The tem¬ 
perature dependence of the isomer shift^ which is shown 
in Fig. 111! is a sum of the temperature-independent 
chemical shift 5c and a temperature dependent contri¬ 


bution 5 r(T) due to the second-order Doppler shift 
5{T) = 5g + 5k{T), (15) 

(16) 

where Mgff denotes the effective mass of the ®^Fe atom. 
i5c can be calculated using i5c = i5(0) — (5r(0). To study 
the influence of the magneto-structural phase transition 
on the isomer shift, the temperature dependence of the 
isomer shift was analyzed in the following way: 5(T^ 
was analyzed in the paramagnetic state with a fixed 0D, 
obtained from Eq. (HI- Then, we extracted the temper¬ 
ature dependence from ^^'^(T'). In a third step we have 
checked whether (5^“(T') can describe the temperature 
dependence of the isomer shift in the magnetically or¬ 
dered state or whether systematic deviations from the be¬ 
havior in the paramagnetic state occur. S^^{T) reveals 
similar behavior for both samples with x = 0.35 and 0.50 
leading to 0.61(3) mm/s and 5™= 0.60(1) mm/s. 
For X = 0.35, below the onset of the magnetic order, 
a deviation from the paramagnetic behavior is found. 
(5™ is reduced to a value of 0.57(1) mm/s. For 

X = 0.50, a reduction to a value of 0.56(1) mm/s 

was observed. Also, this reduction occurs over a wider 


Sr{T) = - 


9 fcg 
16 Meff c 


9d + 8T 
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temperature (K) 


FIG. 11. Temperature dependence of the isomer shift. The 
dashed (solid) lines are a fit in the paramagnetic (magnetically 
ordered) temperature regime using Eq. (1151) with 9d ~ 200 K, 
203 K ,and 213 K for x = 0.35, 0.50, and 0.67, respectively, 
leading to 5^^(r). The deviation from 5^’^(r) corresponds 
to the magneto-structural transition, which causes a change in 
the electron-density at the nucleus. Data and fit for x = 0.67 
are shifted by 5 = — 0.05 mm/s for clarity. 


temperature range as the magnetic transition is broader. 
These reductions in the chemical shifts correspond to an 
increase in the electron density at the nucleus The ori¬ 
gin of this increase can be the structural or the magnetic 
phase transition. The change from a tetragonal to an or¬ 
thorhombic structure changes the lattice parameters and 
hence the volume of the unit cell, which may change the 
chemical shift An increase of the volume and hence 
an increase of the Fe-As distance results in an decrease 
of the electron density at the nucleus corresponding to 
an increased 6c ^ This behavior is for example observed 
in nonmagnetic FeSe, which shows an increase in Sc of 
0.006(1) mm/s at the tetragonal-to-orthorhombic phase 
transition.— A magnetic phase transition may change 
the chemical shift, as observed in metallic iron at the 
Curie temperature, where no structural phase transition 
occurs— The reduction of Jc in metallic iron due to 
the magnetic phase transition is ~ 0.3 mm/s.— This 
may indicate that the origin of the change of the chem¬ 
ical shifts is of both magnetic and structural nature. 
However, previous Mossbauer measurements on the un¬ 
doped compound show contradicting results with either 
no change^ an increase^, or a decrease^ in the chemi¬ 
cal shift. 

For X = 0.67, the temperature dependence of S{T) 
can be properly described by Eq. (Ha and a value of 
dc= 0.56(1) mm/s was obtained. Therefore, no signs of 
a magnetic or structural phase transition were found. 


E. Optimal doped Cao.33Nao.67Fe2As2 

For Cao, 33 Nao. 67 Fe 2 As 2 {x = 0.67), susceptibility 
measurements evidence bulk superconductivity below 


Tc = 34 K, as shown in Fig. [2] Mossbauer spectra down 
to 5 K are shown in Fig. [51 Above 151 K, a doublet struc¬ 
ture with Vzz= 10.3(2) V/A^ is observed. This is con¬ 
sistent with a pure paramagnetic phase. Vzz is constant 
within error bars down to lowest temperatures indicating 
the absence of a structural phase transition. Below 60 K, 
a broadening of the spectra is observed. This broad¬ 
ening indicates small magnetic fields, which were mod¬ 
eled using a Gaussian distribution with a first moment 
(B) = 0. The standard deviation of this Gaussian distri¬ 
bution is constant within error bars, a{B) = 2.2(1) T, 
above Tc and decreases in the superconducting phase 
to 1.7(1) T at 4.2 K. This indicates a competition be¬ 
tween magnetism and superconductivity. Additionally, 
the spectra are nearly symmetric below 60 K and an an¬ 
gle 6 = 55(5)° between the principal axis of the EFG 
and the magnetic hyperfine held is obtained, which is 
close to the magic angle. ZF-/rSR-experiments down to 
5 K are shown in Fig. [3] The time evolution of the muon 
spin polarization exhibits a Gauss-Kubo-Toyabe depo¬ 
larization above 75 K excluding any electronic magnetic 
order. Below 40 K, a weak exponential relaxation sup¬ 
ports short-range magnetic order in a small volume frac¬ 
tion with a MVF smaller than 20 %. By combining both 
local probes, the onset temperature of the weak mag¬ 
netic order is estimated to be 60 K < < 75 K. 

However, room-temperature Mossbauer measurements as 
well as the sharp superconducting transition observed by 
magnetic-susceptibility measurements indicate a homo¬ 
geneous sample. This indicates that the weak magnetism 
is diluted and disordered and persistent even in the op¬ 
timal doping regime, similar to other iron pnictides— “— 

For an investigation of the superconducting phase, TF- 
/iSR measurements were performed in an external mag¬ 
netic helds of fioHext = 11-8 mT perpendicular and par¬ 
allel to the crystallographic c-axis. The magnetic held 
was applied at T>Tc and the corresponding muon spin 
polarization is shown in Fig. I12f a). The weak relaxation 
above Tc is caused by the dipole-dipole interaction of 
the muon spin with randomly distributed dense nuclear 
moments. Additional damping is found in the case of 
a type-B superconductor for hqHci < ^oHext < fJ-oHc 2 
due to the vortex lattice formation. The effect of the 
vortex lattice on the muon spin polarization is shown 
in Fig. [T^ b)-(d). The superconducting signal fraction 
is fully damped after a few fis and « 18 % residual sig¬ 
nal fraction is still oscillating with a precession frequency 
equal to the applied held at times t > 3 /iS. Identifying 
this RslS % signal fraction with the MVF obtained by 
ZF-/rSR measurements, which is of equal value, shows 
that the internal magnetic helds are small compared to 
the 11.8 mT applied held. 

The London penetration depth can be obtained by 
measuring the magnetic-held distribution within the vor¬ 
tex lattice and employing Eq. O- Using the measure¬ 
ments with fioHext\\c, the in-plane penetration depth 
AabCan be directly calculated. For /ioBextTc, contri¬ 
butions from Aac and Abe are measured resulting in an 
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1234 1234 1234 1234 

time (|4S) time (|as) time (|4S) time (|4S) 

a) 38.2 K b) 29.2 K c) 14.1 K d) 1.6 K 


FIG. 12. TF-/iSR spectra for x = 0.67 with ^ollext = 11.8 mT||a6 (upper row) and /ro7/ext=11.8 mT||c (lower row) for temper¬ 
atures above and below the superconducting transition temperature Tc ~ 34 K. The small (Gaussian) damping in spectrum 
(a) is attributed to the dipolar interaction of the muon spin with randomly distributed nuclear moments. The additional 
damping in spectra (b)-(d) is caused by the formation of the vortex lattice in the superconducting state and the associated 
internal magnetic field distribution n{B). It is clearly visible that the damping of the muon precession is stronger in the case 
fJ-oHe^t = 11.8 mT||c. 


effective magnetic penetration depth Ac^fj. Under the as¬ 
sumption of Aq ~ Af,, a value for Ac can be estimated 

• 64 

using** 


Aah — V AflAb Ri Aa —f Ac^fj — V AqAc Ac — , • 

(17) 

The resulting temperature dependence of the in¬ 
verse squared London penetration depth, l/A^f,(T) and 
l/Ac(T), is shown in Fig. fl^ a) and Fig. [T^b), respec¬ 
tively, together with the average internal field (i3), which 
shows a reduction due to the diamagnetic shielding be¬ 
low the superconducting transition. 1/A^(T) was mod¬ 
eled using the phenomenological a-model including two 
independent superconducting gaps with s-wave symme¬ 
try and a fixed Tc = 34 Ki^ The results are shown 
in Tab. [TTl as well as results of ARPES^ {x = 0.67) 
and specific-heat measurements^ (x = 0.68) on single 
crystals. Disorder in the vortex lattice would artifi¬ 
cially reduce the magnetic penetration depth due to the 
broadening of n{B). Therefore, Aab and Ac strictly de¬ 
scribe lower limits (and A“^ an upper limit). Therefore, 



Ai(0)/meV A2(0)/meV lU(Ai) 

A(0)/nm 

Ac 

0.57(8) 

6.7(1.3) 

0.49(4) 

280(46) 

^ah 

0.8(3) 

6(1) 

0.46(8) 

194(17) 

ARPESSi 

2.3 

7.8 



spec. heat^iSii 

2.35 

7.5 

0.75^ 

210(10)ai 


TABLE II. Values of the superconducting gap and the mag¬ 
netic penetration depth obtained by a phenomenological a- 
model analysis of 1/)^{T) and 1/A^(,(T). Ai(0) denotes 

the zero temperature values of the gaps. IT(Ai) and 
IT(A2)=1—VF(Ai) are the corresponding weighting factors. 
A(0) denotes the zero-temperature penetration depth. 


the obtained values for the magnetic penetration depths 
are reduced compared to the values obtained by, e.g., 
specific-heat measurements. To illustrate this effect, the 
temperature dependence of A^^^ heat(^) the cor- 
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(a) 



(b) 



FIG. 13. Temperature dependence of the magnetic penetration depth (a) and (b) (proportional to the superfluid 

density) after field cooling in noH = 11.8 mT, including the ht (solid lines) with a phenomenological a-modelf^ and the average 
magnetic field determined by TF-/rSR. The dashed curve in (b) displays the temperature dependence of \~^ with gap values of 
Ai = 2.35 meV, A 2 = 7.5 meV, W(Ai) = 0.75, and A(0) = 210 nm, obtained by specihc-heat measurements^!^. The reduction 
of the average magnetic field is caused by the diamagnetic shielding of the superconducting phase. 


responding parameter (Ai = 2.35 meV, A 2 = 7.5 meV, 
W(Ai) = 0.75 and A(0) = 210 nm) obtained by specific- 
heat measurements^iffi is plotted in Fig. ll3bl It is clearly 
visible that A’^c. ,heat('^) < for ^ 

This underestimation of the magnetic penetration depth 
in the /iSR experiments may result in a different temper¬ 
ature dependence of A“^(T) and therefore in different gap 
sizes and weighting factors. Additionally, Johnston et al. 
considered the interband coupling of the superconduct¬ 
ing bands and found an intermediate coupling strengthi^ 
The a-model used in this work considers two noninter¬ 
acting superconducting bands, which may also explain 
the different parameter values. Nevertheless, taking into 
account the limitation of magnetic penetration depth 
measurements by means of /iSR experiments and of the 
used a-model, the obtained parameter for Ai(0), A 2 ( 0 ), 
W{Ai), and A(0) are in good agreement with the values 
obtained by ARPES and specific heat experiments 

The anisotropy of the magnetic penetration can be cal¬ 
culated by under the assumption that Ao ~ Ab by^ 


A temperature-independent value of yy = 1.5(4) is the 
smallest observed among the 122 pnictides indicating a 
more 3D-behaviorj^ This behavior is consistent with the 
temperature-independent value of 7 = 1.85(5) for the 
anisotropy of the upper critical fields<^ 


IV. SUMMARY AND CONCLUSIONS 

In summary, we performed muon spin relaxation and 
Mossbauer experiments on Cai_a;Naa;Fe 2 As 2 single crys¬ 
tals with X = 0.00, 0.35, 0.50, and 0.67 resulting in an 
updated phase diagram, which is shown in Fig. 1141 The 
substitution of Ca by Na reduces the onset of the mag¬ 


netic ordering from T™™* = 167(2) K to 161(2) K and 
125(3) K while the magnetic phase transition tempera¬ 
ture width AT = T™®®* — increases from 4 K 

to 21 K and 45 K for x = 0.00, 0.35, and 0.50, respec¬ 
tively. The muon spin precession frequency iz as well 
as the magnetic hyperfine held Rhf, which are propor¬ 
tional to the magnetic order parameter, are reduced as a 
function of the Na-substitution level. Both /rSR as well 
as Mossbauer spectroscopy indicate an increased tilting 
of the magnetic structure upon doping. The magnetic 
phase transition is accompanied by a structural phase 
transition. The lattice dynamics does not change at the 
magneto-structural phase transition. Magnetic suscepti¬ 
bility measurements indicate superconductivity in parts 
of the sample volume for x = 0.35, whereas the sample 
with X = 0.50 shows bulk superconductivity. Therefore, 
as 100 % of the sample is magnetically ordered, coexis¬ 
tence of magnetic order and superconductivity in parts 
{x = 0.35) or in the whole sample (x = 0.50) was ob¬ 
served. A strong reduction of the magnetic order pa¬ 
rameter, as found in the Bai_a;Naa;Fe 2 As 2 series^ is 
not observed for Cai_xNaa;Fe 2 As 2 with x = 0.35. For 
X = 0.50, a small reduction of ~ 7 % was observed. We 
applied a Landau theory to describe the reduction of the 
magnetic order parameter showing that the magnitude of 
the reduction depends on the coupling strength and the 
Tc/Tjq ratio. A linear relation between the reduction of 
\Mco\'^/\M q\’^{T = 0 ) on Tc/Tat has been found for several 
superconducting 122 -pnictide systems that microscopic 
coexistence. 


For X = 0.67, diluted and weak magnetism below 60- 
75 K as well as bulk superconductivity with Tc = 34 K is 
found. The s-wave symmetry of the two superconducting 
gaps as well as the value of the larger gap agrees well with 
recent ARPES-and specific-heat measurementsj ^^i^^ 
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